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Measurement of Parity Violation in the Early Universe using Gravitational-wave 

Detectors 
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A stochastic gravitational-wave background (SGWB) is expected to arise from the superposition of 
many independent and unresolved gravitational-wave signals, of either cosmological or astrophysical 
origin. Some cosmological models (characterized, for instance, by a pseudo-scalar infiaton, or by 
some modification of gravity) break parity, leading to a polarized SGWB. We present a new technique 
to measure this parity violation, which we then apply to the recent results from LIGO to produce the 
first upper limit on parity violation in the SGWB, assuming a generic power-law SGWB spectrum 
across the LIGO sensitive frequency region. We also estimate sensitivity to parity violation of the 
future generations of gravitational-wave detectors, both for a power-law spectrum and for a model 
of axion inflation. This technique offers a new way of differentiating between the cosmological and 
astrophysical sources of the isotropic SGWB, as astrophysical sources are not expected to produce 
a polarized SGWB. 

PACS numbers: 



Introduction. — A stochastic gravitational-wave back- 
ground (SGWB) is expected to arise from the superpo- 
sition of gravitational waves (GWs) from many uncor- 
rected and unresolved sources. Numerous cosmological 
SGWB models have been proposed, including inflation- 
ary models models based on cosmic (super)strings 
[El, @ i and models of alternative cosmologies Q ■ Further- 
more, various astrophysical models have been proposed 
based on integrating contributions from astrophysical ob- 
jects across the universe, such as compact binary coales- 
cences of binar y n eutron stars and/or black holes [1, 



magnetars 



ii 



or rotating neutron stars [1 
eral searches for the unpolarized isotropic II 31-41 




anisotropic SGWB [T7J have been conducted using 



data acquired by interferometric GW detectors LIGO 
3 1|| and Virgo [2(J. These searches have established 
upper limits on the energy density in the SGWB, and 
have started to constrain some of the proposed models 

In this Letter we present the first upper limits on the 
circularly polarized isotropic SGWB. Since astrophysi- 
cal sources are unlikely to produce a circularly polar- 
ized isotropic SGWB, detecting polarization asymmetry 
in the SGWB is potentially an excellent way of distin- 
guishing the cosmological component from the possibly 
dominant astrophysical one. Such asymmetry could be 
generated, for instance, through the gravitational Chern- 
Simons term [22I [23j , from the imaginary part of the Im- 



mirzi parameter [24j |. in some power-counting renormal- 
izable theories of gravity [251 ] . or, as we discuss in more 
detail, from an axion infiaton [26j. We follow the for- 
malism developed in |27| . modified to address a polarized 
SGWB as discussed in [23| . Using the latest SGWB mea- 
surement with LIGO detectors [15| , we apply this formal- 
ism to produce the first constraints on parity violation 
for a generic power-law SGWB spectrum. We also esti- 
mate the sensitivity of the upcoming second-generation 
GW detectors to the power-law and axion-inflation par- 
ity violating models — Advanced LIGO (aLIGO) [IH de- 
tectors at Hanford, WA (HI) and Livingston, LA (LI), 
Advanced Virgo [H in Italy (VI), GEO-HF in Ger- 
many, and KAGRA [3l|,[32| m Japan (Kl) are expected 
to have ~ 10 x better strain sensitivities than the first- 
generation detectors, and to produce first science-quality 
data in 2015. Finally, we consider an example configu- 
ration of a pair of third-generation GW detectors, with 
strain sensitivity similar to the proposed Einstein Tele- 
scope 33 1 . 

Search Formalism. — We start from the plane-wave ex- 



pansion of the metric at time t and position x [231 127 



(t,x) = J2 / df dClh A (f,Cl)e 



-2-nif{t—x-iX) A 



(1) 

where e^ b (f2) is the polarization tensor associated with a 
wave traveling in the direction f2, and / is frequency (we 
use natural units c = h — 1). We consider the left- and 
right-handed correlators [23j |: 
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where hi, = (h + + ih x )/y/2, hR = (h + — ih x )/V2, and 
+ and x arc the standard plus and cross polarizations. 

Note this is the point of departure from the past 
searches for unpolarizcd isotropic SGWB, which assume 
V = 0. Further note that (hnh* L } vanishes due to statis- 
tical iso tropy. We then compute the normalized energy 
density daS: 



now(/) 



/ dp 



GW 
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1(f) 
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where dpow is the energy density in the range [/, f + df], 
Gjv is Newton's constant, and p c the critical energy den- 
sity of the universe Q. We also compute the standard 
cross-correlation estimator [27} : 

/+oo r + co 

df / df'6 T (f - f')((sl(f)s 2 (f'))Q(f) 
-oo J — oo 

3H$T f d/ Q GW (/)7z(/)Q(/) ? (4) 



where 



^Gw(/)7/(/) = Ogw(/) [7/(/) + Whv(f)} (5) 
7/ (/) = JL. J dCl(F+F+* + F?*)e 2niffl - Ai 

7v (/) = - A /" d(li{F+F** - F*F+*)e 27Ilih - AS . 



Here, T is the measurement time, <5t(/) = 
sin(7r/T)/(7r/), si(/) and S2(/) are Fourier trans- 
forms of the strain time-series of two GW detectors, 



d„ is the contraction 



Q(f) is a filter, and F 7 f 
of the tensor mode of polarization A, e^ b , with the 



J ab n 



response of the detector n, d£ 6 I45j. The factor 7/(/) 
is the standard overlap reduction function arising from 
different locations and orientations of the two detectors, 
and jv{f) is a new function, associated with the parity 
violating term. Figure [1] shows these functions for two 
real detector pairs. Finally, IT(/) = V(f)/I(f) encodes 
the parity violation, with maximal values II = ±1 cor- 
responding to fully right- or left-handed polarizations. 
Setting IT = reproduces the standard unpolarized 
SGWB search ^ . 

Assuming stationary Gaussian detector noise (uncorre- 
lated between two detectors), the estimator for the vari- 
ance associated with Y is 12711: 



^ = 1 / 



(0) 



where P n (f) are the one-sided noise power spectral den- 
sities of the two GW detectors. In practice, we divide 
the sensitive frequency band of the GW detectors into 
bins A/ = 0.25 Hz wide [15J. We then compute the es- 
timator Yi and the variance af for each frequency bin 
i assuming a frequency independent spectrum template 
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FIG. 1: Overlap reduction functions for detector pairs Hl-Ll 
and Ll-Kl. 



S]gw(/) = for each bin. Optimization of the signal- 
to-noise ratio then leads to the following optimal filter for 
a frequency-independent GW spectrum in the frequency 
bin ft [13: 



Q(fi) = 



liifi) 



ffPxifiWfi. 



(7) 



with normalization constant M chosen so that (Yi) = Q.q. 

To perform parameter estimation in the parity violat- 
ing models, we adopt the Bayesian approach introduced 
in 21] . We define the following likelihood function: 



L(Yl, <Ji\9) oc exp 



2 ^ 



(Y-n> M (f i; 6))< 



(8) 



where the sum runs over frequency bins fi, and Q,' M (/$; 9) 
is the SGWB energy density spectrum specified by the 
free parameters 8, multiplied by the parity-violating cor- 
rection (l + n(/)7y(/)/ 7/ (/)), as defined in Eq. [5] Mul- 
tiplying the likelihood with the prior distribution for 9 
yields the Bayesian posterior distribution, which can then 
be used to extract confidence intervals for 9. In the sub- 
sequent results, we take all priors to be flat within the 
plotted range of 9 (and zero elsewhere). 

Power Law. — We apply this formalism to the power- 
law spectrum, 51m(/) = ^a(///re/) a , which is appropri- 
ate for most SGWB models in the LIGO frequency band. 
We set f re f — 100 Hz, and assume n(/) = n = const. 

The results are shown in Figure [5] The top-left plot 
of the Figure shows the 95% CL exclusion curves in the 
n — Q a plane for several values of the spectral index a, 
using the latest SGWB measurement with LIGO detec- 
tors (specifically, using the data plotted in Figure 6 of 
fl5l| ) . This plot represents the first constraint on parity 
violation in the early universe using GW observations. 
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FIG. 2: Top-left: 95% CL limit curves based on the most recent LIGO result [l5| . for several values of a. Top-right: Expected 
sensitivities (at 2a level) for the advanced LIGO Hl-Ll pair and for an example of a third-generation detector pair (see text 
for more detail). Bottom-left: 95% CL contours are shown for several examples of simulated SGWB with parity violation, 
assuming a — 0, standard strain sensitivities, and 1 year of observation. The x's denote the signal simulation parameter values. 
The lightest-gray line corresponds to the recovery with Hl-Ll, medium-gray line to the recovery with H1-L1-V1, and the black 
line to the recovery with H1-L1-V1-K1 network. Bottom-right: assuming a SGWB with maximal parity violation (II = +1), the 
lines denote Q a needed for a given a to detect the SGWB and exclude LI = at 95% confidence, using two second-generation 
detector networks and the example of a third-generation detector pair. 



The top-right plot of Figure [5] shows the expected sensi- 
tivities for the detection of a SGWB at a 2a level of the 
Advanced LIGO Hl-Ll detector pair, assuming standard 
expected strain sensitivity (28[. It also shows the sensi- 
tivity of an example configuration of two third-generation 
detectors, assuming two L-shapc detectors of strain sen- 
sitivity similar to the Einstein Telescope design j3^|, lo- 
cated at opposite poles of the Earth, with respective arm 
orientation rotated by n/8. Such a pair of detectors has 
not been proposed yet, so this example is only meant to 
illustrate the potential reach of third-generation detec- 
tors, in a configuration that suppresses neither 77 nor 
jy terms. In both cases we assume 1 year of observa- 
tion. As expected, the second-generation detectors will 
provide 3-4 orders of magnitude better SGWB measure- 



ment as compared to the first-generation detectors, and 
the third-generation detectors could provide yet another 
3-4 orders of magnitude improvement. 

The bottom- left plot of Figurc[2]illustratcs the recovery 
of several simulated parity violating SGWB signals using 
different networks of second-generation GW detectors. 
Having multiple detector pairs in the network does not 
significantly improve the theoretical uncertainty <Ji (Eq. 
[5]) over the Hl-Ll detector pair. However, multiple de- 
tector pairs introduce multiple 7y's in the search, which 
helps break the degeneracy between the parity conserv- 
ing and violating terms in Eq. [5J As a result, the 95% CL 
contours are significantly tighter in n when using multi- 
ple detector pairs. Finally, the bottom-right plot shows 
the amplitude VL a required to detect parity violation 
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with the Hl-Ll and H1-L1-V1-K1 second-generation de- 
tector networks, and with the third-generation example 
pair (assuming maximal parity violation in the SGWB, 
n = +1). In particular, for a = the H1-L1-V1-K1 
detector network could detect parity violation at ampli- 
tudes 4x smaller than the Hl-Ll pair alone, while the 
third-generation detectors could detect it at amplitudes 
another ~ 4 orders of magnitude smaller. 

Parity violation during inflation. — We now apply the 
above formalism to a specific example of a parity violat- 
ing SGWB. We concentrate on an example characterized 
by standard gravity. Specifically, we consider a natu- 
ral vast class of inflationary models where the flatness 
of the potential V of the inflaton <fi is protected by an 
approximate shift symmetry <p ~ > 4> + constant. The 
most studied realization of this is the axion field 34j . 
and many models, starting from [35j . have employed this 
idea in the context of inflation. Generally, after infla- 
tion an axion inflaton predominantly decays into gauge 
quanta through the coupling AC ~ j^cpFF, where f$ is 
the axion decay constant, C is a parameter which is nat- 
urally of order one, F is the field strength of the gauge 
field (for simplicity, a U(l) field is considered here), and 
F its dual. This interaction leads to an interesting gauge 
field production also during inflation [36j . as it modifies 
the dispersion relation of the two gauge field helicities to 



C(j) 



(9) 



Here, a is the scale factor of the universe that grows 
nearly exponentially during inflation, a ~ e , while the 
other dynamical quantities are slowly varying. Any gauge 
field mode is identified by its constant comoving momen- 
tum k. The second term in ui± becomes dominant as the 
Compton wavelength A ~ a/k becomes greater than the 
Hubble horizon H^ 1 (this moment is dubbed "horizon 
crossing"; we denote by the value of £ at this mo- 
ment). As a consequence, one of the two helicities expe- 
riences a tachyonic growth. For definitcness, we assume 
£ > 0, so that the growing helicity is the + one. Modes of 
progressively smaller A (larger k) cross the horizon later 

during inflation. Since £ cx grows during inflation, £& 
is a slowly increasing function of k. Since the amplitude 
of the produced gauge modes is cx e 71 ^ J3(|, even this 
slow increase can result in a much stronger production 
at short wavelengths. 

The produced gauge quanta generate inflaton pertur- 
bations and GW through A + A + — > 8<j> and ^4+A + — > h 
processes These "sourced" signals superimpose in- 
coherently to the standard "vacuum" inflaton and met- 
ric perturbations. For example, the energy density in 
GW generated during inflation is flow oc (Vacuum) + 
(^-sourced); an d identically for 5<fi. The A + ^4 + — ► h in- 
teraction mostly produces GW of left chirality [2(| , and 
therefore n ~ — 1 whenever /i snllrm d 3> h vacumn (equiv- 



Detector Network 


P 


£ (sensitivity) 


£ (exclude II = 0) 


2 nd Gen H1-L1-V1-K1 


1 


2.3 


> 2.66 


2 nd Gen H1-L1-V1-K1 


2 


2.2 


2.6 


3 rd Gen 


1 


1.8 


2.0 


3 rd Gen 


2 


1.9 


2.0 



TABLE I: Column 3 shows the 2<r sensitivity to £ for different 
future GW networks (column 1) in the axion-based inflation 
model discussed in the text, for two power-law inflaton po- 
tentials: p = 1 or p — 2 (column 2). Column 4 indicates the 
value of £ needed in order for the respective detector network 
to detect parity violation (reject II = 0) at 2a. 



alently, n ~ +1 for £ < 0). The 5<^ SO urced signal is 
highly non-gaussian [33, [38| , while the primordial density 
perturbations are highly gaussian. Therefore, S(f> sourcc d 
needs to be subdominant to the vacuum signal, which, 
according to the estimate of [HI , results in the constraint 
£cmb 2.66 at the moment in which the modes of cos- 
mological scales (those that determine the Cosmic Mi- 
crowave Background anisotropics) were produced flrjj ]. 



< h 



vacuum 



When this limit is respected, h so ^ 

mological scales, and is unobservable 26[ 37, 38 in the 
minimal implementations of the mechanism. However, 
given the increase of £& with k, the produced A + and, 
consequently, the amount of h soulcc d steeply increases 
with increasing frequency [39j , to a level that can be ob- 
served by terrestrial GW detectors. The computation of 
this effect requires taking into account the backreaction 
of the produced gauge fields on the late inflaton dynam- 
ics |40j |. and leads to the results summarized in Figures 
5 and 6 of that work. It was shown there that the signal 
can be observed by Advanced LIGO for values of £ that 
do not conflict with the CMB bound. 

We apply the formalism developed above to this in- 
flation model. In particular, we consider two power-law 
inflaton potentials, V ~ <fi p for p = 1 (typical of axion 
monodromy |4lj |) and p = 2 (a generic Taylor expan- 
sion around the minimum of a potential). For each p, 
^m(/;£) in Eq. [8] is a function of the single parameter 
£. We then compute the expected 2a sensitivities to £ 
for the second-generation H1-L1-V1-K1 detector network 
and for the example third-generation detector pair, as- 
suming flat prior distribution for £ up to the CMB bound 
of [38| . The results are shown in Table |TJ The second- 
generation detectors will begin to constrain the £ pa- 
rameter beyond the CMB gaussianity constraint of j38| . 
while the third-generation detectors may reach £ ~ 1.8. 
Somewhat larger values of £ are required in order to de- 
tect parity violation in the observed SGWB and rule out 
n = at the 2cr level. 

We note that £ can also be constrained by taking ad- 
vantage of its impact on J^sourccd, which leads to the 
non-gaussian signature mentioned above, to the growth 
of the power spectrum of the scalar perturbations at the 
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highest values of k that can be observed in the CMB |42J , 
and to the production of primordial black holes |43| . The 
bounds from the first two effects [42j are weaker than 
those from the SGWB. Those from the primordial black 
holes are potentially stronger, but they carry an uncer- 
tainty associated with how many scalar perturbations are 
produced at such small scales, and with how many black 
holes they create 43|. No such uncertainties are present 
for the SGWB signal that we have studied here. 

Conclusions. — We introduce a new formalism for mea- 
suring polarization asymmetry in the SGWB, which both 
improves the sensitivity of the SGWB searches and of- 
fers a new handle for distinguishing between the cos- 
mological and astrophysical SGWB sources. We apply 
the formalism to the recent LIGO data to produce the 
first constraints on the parity violation in the SGWB. 
We also estimate the sensitivities of future GW detectors 
to the parity violating SGWB signals in the framework 
of a generic power-law SGWB spectrum and in the spe- 
cific example of an axion-based inflation model, demon- 
strating the substantial advantage of networks allowing 
multiple non-collocated detector pairs. 
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